Electron transfer theory revisit: Quantum solvation effect 
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The effect of solvation on ET rate processes is investigated on the basis of the exact theory 
constructed in J. Phys. Chem. B 110, xxx (2006). The nature of solvation is studied in a close 
relation with the mechanism of ET processes. The resulting Kramers' turnover and Marcus' inversion 
characteristics are analyzed accordingly. The classical picture of solvation is found to be invalid when 
the solvent longitudinal relaxation time is short compared with the inverse temperature. 



I. INTRODUCTION 

Chemical reaction in condensed phases is intimately 
related to the Brownian motion in solution. Einstein's 
paper on Brownian motion! showed the first time the 
fluctuation-dissipation relation (FDR). The fluctuations 
of surrounding molecules are responsible for both ag- 
itation and friction on the Brownian particle. These 
stochastic events of energy exchange between system and 
bath lead eventually to thermal equilibrium. Brownian 
motion is characterized by the stochastic force. The FDR 
leads to a Fokker-Planck equation, which is the classical 
reduced equation of motion that governs the Brownian 
motion at an ensemble average level. This approach has 
been exploited by Kramers in his construction of iso- 
merization reaction rate theoryi£ The resulting rate is 
shown to have a maximum in an intermediate viscos- 
ity region. This celebrated Kramers' turnover behavior 
clearly demonstrates the dual role of solvent on reaction 
rateiSi 3 . 

Electron transfer (ET) is the simplest but a 
pivotally important chemical reaction system. It 
constitutes another class of systems whose depen- 
dence on solvent environment has been extensively 
studied since the pioneering work by Marcus in 
1950s, 4 ! 5 ! 6 ! 7 ! 8 ! 9 . 10 ! 11 . 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 ! 19 ! 20 However, it is 
often treated in different way from the traditional chem- 
ical reaction involving bond breaking and/or formation. 
In the latter case, either the equation of motion for a 
particle over the barrier or the flux-flux correlation func- 
tion approach on the basis of the transition-state the- 
ory is usedi 3 i 21 ' 22 i 23 The standard treatment in the ET 
research field is rather a type of transfer coupling cor- 
relation function formalism, based on the assumption 
that the nonadiabatic coupling matrix element V is not 
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Depicted in Fig.^ is the schematics of an elementary 
donor-acceptor ET system. Here, E° denotes the re- 
action endothermicity; V a (Vf>) represents the potential 
surface of the solvent environment for the electron in the 
donor (acceptor) state; U = Vj, — V a — E° is the solvation 
coordinate; while A = (U) is the solvation energy, with 
(■ • •) denoting the initial bath ensemble average. At the 



crossing (U + E° = 0) point, V a = V b = (E° + A) 2 /(4A) 
that amounts to the ET reaction barrier height. The 
celebrated Marcus' rate theory read o 4 i 5 i 6 



k = 



V 2 /h 



exp 



(E° + A) 2 



4Afc B T 



(1) 



It is a classical Franck-Condon theory, assuming that the 
solvent relaxation is much slow compared with the elec- 
tronic transition. Exploited in Eq. Q is also the classical 
FDR: {U 2 ) - {U) 2 = 2k B T{U). Quantum extension of 
Marcus' theory has been formulated in the weak trans- 
fer coupling regimens . The dynamic solvation effect is 
introduced by the solvation correlation function, 



C(t) = (U(t)U(0)) - (Uf 



(2) 



Nonperturbative rates have also been formulated on 
the basis of fourth-order transfer correlation func- 
tions, followed by certain resummation schemes to par- 
tially account for the nonperturbative transfer coupling 
effects, 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 The resulting rates, despite 
of the resummation approximation involved, do recover 
the celebrated Kramers' turnover behavior. 

The main purpose of this work is to elucidate some 
distinct solvation effects on the ET rate processes. The 
quantum nature of solvation arises from the fact that the 
solvation coordinate is an operator and its correlation 
function must be complex. The elementary ET system in 
Debye solvents will be studied on the basis of the exact 
and analytical rate theory developed in Ref. 24, which 
will be referred as Paper I hereafter. Section [H] sum- 
marizes the theoretical results of Paper I. This is a re- 
duced quantum equation of motion based formalism; i.e., 
the quantum version of Kramers' Fokker-Planck equation 
approach. The key quantity now is the reduced density 
matrix, p(t) — tr B pT, defined as the trace of the total 
system-bath density matrix over the bath subspace. Nu- 
merical results will be presented and discussed in Sec. lIIII 
Finally, Sec. II VI concludes this paper. 
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II. AN EXACT AND ANALYTICAL THEORY 

This section summarizes the exact and analytical rate 
theory, developed in Paper I, for the ET in Debye solvents 
at finite temperatures. Let us start with the following 
form of the reduced Liouville equation, 



p(t) = --[H, p(t)} -j dTll(t, r)p(r). (3) 

For the present ET system (Fig.^), the reduced system 
Hamiltonian reads 



H=(E°+ X)\b)(b\ + V(\a)(b\ + \b)(a\). 



(4) 



This is a time-independent system, for which the dissi- 
pation memory kernel LT(t, r) = Il(t — t). As a result, 
Eq. (0 can be resolved in its Laplace domain as 



sp(s) - p(0) = -iLp(s) - n(s)p(s). 



(5) 



Here L = 7i _1 [.H~, •] is the reduced system Liouvillian. 

A simplification arises for the ET system in Debye sol- 
vents at finite temperature. The solvation correlation 
function assumes now (for t > 0) 



C(t) = A(2fc B T - ihj)e 



-it 



(6) 



Here, tl = 7 _1 denotes the longitudinal relaxation time 
of the Debye solvent. In this case, Eq. (0 can be for- 
mulated exactly in terms of a continued fraction Green's 
function formalism. Let H(s) = (s) and 



s + ih + n(™) (s 



n > 0. 



(7) 



The continued fraction hierarchy is now the relation be- 
tween n ( ")(s) and g(™ +1 )(s + 7); cf. the eq (16) of Pa- 
per I. For the elementary ET system subject to the De- 
bye longitudinal relaxation, it is found that U\ n >, which 
is Hermite satisfying n^™, fefe/ = n^"^, fc , has only three 
nonzero tensor elements together with their complex con- 
jugates. Denote the three nonzero tensor elements of n> n ) 
as 



X = UbaMai V = nfca.ab, Z = UbaMb- 



(8) 



Implied here, and also whenever applicable hereafter [cf. 
Eqs. (|10afl - ljllbjl ]. are the superscript index (n) and ar- 
gument s, if they are the same in the both sides of the 
individual equation; otherwise they will be specified. The 
continued fraction hierarchy that relates n*-™' (s) with 
g(n+i) ^ s _|_ ^ can now ex p resgec | i n terms of 

x (n) (s) = r 1 (n+ l)X( n+1 )(s + 7 ), (9a) 
= -7y*(n + l)r( n+1 )( S + 7 ), (9b) 
z {n) {s) = {■q-i 1 *){n + l)Z {n+1 \s + 1 ). (9c) 

Here, {X,Y,Z} {n \s) are the counterpart tensor elements 
of the Green's function Q^ n \s). Their relations to the 



nonzero elements {x, y, z]»(s) of IlW(s) via Eq. Q 
can be evaluated analytically on the basis of the Dyson 
equation technique. The final results are 



X — GbaMa — 
y = Gba.ab = 



a* +f3* 



P-v 



(10a) 
(10b) 



Z = GbaM = -~[{z-iV/ti)X + (z* + iV/h)Y](10c 



with 



P = 



s+(i/h)(E° 
s-\V/hf[2- 



- A) + x, 
ihz/V]. 



(11a) 
(lib) 



The kinetics rate equation can be readily obtained via 
Eq. (0 by eliminating the off-diagonal reduced density 
matrix elements. It leads to a linear algebraic equation in 
the Laplace domain that corresponds to the generalized 
rate equation with memory rate kernels in time domain. 
The resulting ET rate resolution reads as [the eq (36a) 
of Paper I] 



k(s) 



2\V\ 2 a(s)+y(s) 

V H*)l 2 - \v{s) 



(12) 



The rate constant k = fc(0) that will be numerically stud- 
ied in the next section amounts to the time integral of 
the memory rate kernel in the aforementioned generalized 
rate equation. 

As analyzed in Paper I, the infinity inverse recursive 
formalism [Eqs. 10 and (|10fl ] can be truncated by set- 
ting {x, y, z}( N } = at a sufficiently large anchoring N. 
The resulting {x,y 7 z}^ = {x,y,z} that are required 
by Eq. ifl^) are exact up to the (27V) th -order in the 
system-bath coupling. The convergence is guaranteed 
also via the mathematical continued fraction structure 
involved. Apparently, if the rates are needed at a speci- 
fied s', one shall start with {x, y, z}^} 3 , +N = 0. The 
backward-recursion relations, Eqs. @ with Eqs. ltTU|l . 
will then lead to {x,y,z}^ s ,^ N _ 1 ^ 1 and so on, until 

{x, y, z}^} s , = {x, y, z} s=s > are reached for evaluating the 
required k(s') [Eq. (|12fl ]. The above reduced dynamics- 
based ET rate formalism is exact for the Debye solvents 
in finite temperatures. However, the FDR, which relates 
the real and imaginary parts of the solvation correlation 
function, is adopted in Eq. © in a semiclassical manner. 
As a consequence the reduced density matrix and rates 
may become negative if the temperature is too low. 



III. QUANTUM SOLVATION EFFECTS: 
NUMERICAL RESULTS 



We are now in the position to elucidate some distinct 
solvation effects on the ET reaction rate k = k(s = 0) [cf. 
Eq. (|12|) ]. Numerical results will be presented in relation 
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to the celebrated Kramers' turnover and Marcus' inver- 
sion behaviors, exemplified with the ET reaction systems 
of V = 1 kJ/mol and A = 3 kJ/mol at T = 298 K. 

It is noticed that the solvation longitudinal relax- 
ation time tl is considered proportional to the solvent 
viscosity^&ii The Kramers' turnover characteristics can 
therefore be demonstrated in terms of the rate k as a 
function of the scaled solvent relaxation time T^/r^ev 
Here, Tthcr = fo/{kB.T) denotes the thermal time, which 
at the room temperature is about 26 fs. In the Debye sol- 
vent model of Eq. JJjJl, the quantum nature of solvation 
enters via the semiclassical FDR that relates the real and 
imaginary parts of the correlation function. In contrast, 
the classical solvation is characterized by the real part 
only. As Im?7/Rc?7 = — 0.5Tther/'?L, it is anticipated that 
the quantum nature of solvation can only be prominent 
in the low viscosity (tl < T t her) regime. It is also con- 
sistent with the physical picture that the high viscosity 
(or slow motion) implies a large effective mass and thus 
leads to the classical solvation limit. 

Figure[5]depicts the rates k as functions of T^/T t hei for 
two typical systems, being of the endothermicity values of 
E° = —A and E° = 0, respectively. Observed in the high 
viscosity (tl /Tther > 1) regime for each of the systems is 
the celebrated Kramers' fall-off behavior.?' 3 This is the 
well established classical solvation picture of the diffu- 
sion limit: the higher the solvent viscosity is, the more 
backscattering (or barrier-recrossing) events will be. The 
fact that kE°=-x > &e°=o observed in the high viscos- 
ity regime is also anticipated from the classical solvation 
picture [cf. Fig^or Eq. QJ]: That E° = —A represents a 
classical barrierless system where the celebrated Marcus' 
inversion takes place. 

In the low viscosity (tl/ T tri cr < 1) regime the clas- 
sical picture of solvation is however invalid. The ob- 
served rate in the symmetric (E° = 0) system, is appar- 
ently tunneling dominated due to Fermi resonance. The 
most striking observation is that the so called barrierless 
(E° — —A) system exhibits now clearly the Kramers' 
viscosity-assisted barrier-crossing characteristics in the 
low viscosity regime. This suggests that there is an effec- 
tive barrier for the ET system with the classical barrier- 
less value of E° + A = 0; this effective barrier is viscosity 
dependent and vanishes as tl increases. 

Now turn to the Marcus' inversion characteristics for 
the rate A; as a function of reaction endothermicity E°. 
Depicted in Fig.[3]are the resulting inversion curves, with 
TL/Tthcr = 0.1,1, and 10 to represent the low (solid- 
curve), intermediate (dot-curve), and high (dash-curve) 
viscosity regimes, respectively. In the classical solvation 
picture the inversion occurs at E° = —A, as it represents 



a classical barrierless system. This picture is only valid 
in the high viscosity regime; see the dashed curve with 

TL /Tthcr = 10. 

In the low viscosity regime, according to the analysis 
presented for Fig.0 there is always a nonzero barrier for 
the ET reaction, covering over the entire range of E° 
including the value of E° = —A. This explains the inver- 
sion behavior of the solid curve in Fig. [3] that is peaked 
only at the resonant position of E° — 0. As the viscosity 
increases, the inversion region smoothly shifts from the 
resonant peak position E° = to the classical barrierless 
position of E° = —A. 

To explain the asymmetric property of the inversion 
behavior observed in Fig.EJ let us recall that k(—E°) 
amounts to the backward reaction rate and k{E°) < 
k(—E°) for an endothcrmic (E° > 0) reaction. This 
leads immediately to the asymmetric property of the solid 
curve in Fig-El i n which the blue (endothermic) wing falls 
off faster than its red (exothermic) wing. The degree of 
asymmetry decreases as the viscosity increases. Only in 
the high viscosity regime does the inversion curve behave 
classically, which is symmetric (but may not be parabolic 
unless the transfer coupling V is small) around its clas- 
sical inversion position of E° = —A. 



IV. SUMMARY 

In summary, we have investigated in detail the effect 
of solvation on ET rate processes. The nature of solva- 
tion is studied in a close relation with the mechanism of 
ET processes in terms of Kramers' turnover and Marcus' 
inversion characteristics. The classical picture of solva- 
tion is found to be invalid in the low viscosity regime, 
which can be well measured by the scaled longitudinal 
relaxation time of T-^/rthcn where Tthcr = f^l{k^T) is the 
thermal time. The present study is carried out on the 
basis of the exact rate theory for the simplest ET system 
with a single solvation relaxation time scale. Neverthe- 
less, the basic results obtained here are expected to be 
valid to a general ET system in a realistic solvent envi- 
ronment of multiple relaxation time scales. 
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FIG. 1: Schematics of solvent environmental potentials V a 
and Vb for the ET system in the donor and acceptor states, 
respectively, as the functions of the solvation coordinate U = 
Vb — V a — E° , with E° being the ET endothermicity and A = 
(U) the solvation energy. The classical barrierless system is 
that of E° + A = 0. 



FIG. 2: ET rates (k) as functions of scaled solvent longitu- 
dinal relaxation time (Tx/Tthor) for the symmetric (E° — 0) 
and the classical barrierless (E° + A = 0) systems at T=298 
K, with A = 3 kj/mol and V = 1 kj/mol. The thermal time 
Tthcr = fr/iksT) = 10~ 1,6 ps at the room temperature. 



FIG. 3: ET rates (k) as functions of reaction endothermicity 
(E°). Three values of the relative relaxation time, TL/Vthcr = 
0.1,1,10 are used to represent low, intermediate, and high 
viscosity regimes. 
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